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ABSTRACT 

We study the polarization properties of relativistic reconfinement shocks with chaotic mag- 
netic fields. Using our hydrodynamical model of their structure, we calculate synthetic polar- 
ization maps, longitudinal polarization profiles and discuss the spatially averaged polarization 
degree as a function of jet half-opening angle Qj, jet Lorentz factor Tj and observer inclina- 
tion angle to the jet axis 8 {, s . We find, that for o b s < Qj the wave electric vectors are parallel 
in the vicinity of the structure ends and perpendicular in between, while for o t, s > Qj the 
polarization can only be perpendicular. The spatially averaged polarization degree does not 
exceed 30%. Parallel average polarization, with polarization degrees lower than 10%, have 
been found for 9 t, s < Qj under the condition TjQj > 1. As earlier works predicted the par- 
allel polarization from relativistic conical shocks, we explain our results by discussing conical 
shocks with divergent upstream flow. 

Key words: polarization - shock waves - galaxies: jets. 



1 INTRODUCTION 

A great variety of polarization properties has been found in rel- 
ativistic AGN jets on different wavelengths (radio, mm and opti- 
cal) and scales (subparsec through kiloparsec). From the theoretical 
point of view, high polarization degrees in this wavelenght range 
implicate the synchrotron emission mechanism taking place in the 
presence of ordered magnetic fields. The polarization electric vec- 
tors are being observed parallel or perpendicular to the jet axis. 
Perpendicular polarization, prevailent in the large-scale jets in ra- 
dio (Bridle et al. 1994) and optical (Perlman et al. 2006) bands, has 
been particularily difficult to explain, as it requires the dominance 
of longitudinal magnetic fields in the emission regions. In expand- 
ing jet the parallel (poloidal) component of the magnetic field de- 
cays faster than the perpendicular (toroidal) one, so one needs a 
process, in which poloidal component is amplified on large scales 
(Begelman et al. 1984). The velocity shear at the jet boundary is 
usually invoked (Laing 1981). 

Chaotic (tangled) magnetic fields are thought to dominate the 
jets at distances larger than a few parsecs, but they must be statis- 
tically anisotropic to produce a net linear polarization (Korchakov 
& Syrovat-Skii 1962). It has been pointed out by Laing (1980), 
that initially isotropic chaotic magnetic field becomes anisotropic 
after crossing the shock front due to compression of plasma. Po- 
larization degree of emission from compressed magnetic field has 
been calculated by Hughes et al. (1985) and more general for- 
mulae, incorporating distortions due to velocity shear, have been 
given by Matthews & Scheuer (1990) and Laing (2002). Polariza- 
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tion from stationary relativistic conical shocks has been studied by 
Cawthorne & Cobb (1990), hereafter CC90. They found that high 
degrees of parallel polarization may be obtained, but the degree of 
perpendicular polarization is limited to ~ 10%. Since the knots of 
blazar jets are observed with a perpendicular polarization of higher 
degree (e.g. Marscher et al. 2002), Cawthorne (2006) introduced 
large-scale poloidal magnetic field component in the upstream flow. 

In weakly magnetized AGN jets, a sequence of so-called re- 
confinement shocks forms, resulting from the interaction between 
the jet and the external medium (Sanders 1983). In our previous 
paper (Nalewajko & Sikora 2009, hereafter NS09) we studied the 
structure and energy dissipation efficiency of axisymmetric recon- 
finement shocks and compared our results to the analytical formu- 
lae given by Komissarov & Falle (1997). We will now use a model 
developed there as a basis for calculating the polarization of emis- 
sion originating behind the shock front from chaotic magnetic fields 
compressed at the shock. A reconfinement shock is treated as a set 
of conical shocks with inclination angle dependent on the position 
along the symmetry axis. The upstream flow is assumed to be ex- 
panding freely, thus it diverges. We begin with demonstrating the 
difference between conical shocks with parallel and divergent up- 
stream flow. 

This paper is organized as follows. In §2 we present the 
scheme for calculating the degree and positional angle of linearly 
polarized emission from oblique shocks. In §3 we study the polar- 
ization from conical shocks with diverging upstream flow. In §4 we 
study the polarization from reconfinement shocks: the polarization 
maps are presented in §4.1, the longitudinal polarization profiles in 
§4.2 and the spatially averaged polarization degrees in §4.3. Our 
results are discussed and summarized in §5. 
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We use a term 'perpendicular/parallel polarization', meaning 
the orientation of the wave electric vectors with respect to the jet 
axis. The primed quantities are those measured in the frame co- 
moving with the downstream plasma, in contrast to the quantities 
measured in the external frame. 



2 CALCULATING THE POLARIZATION FROM THE 
OBLIQUE SHOCK 

In weakly magnetized shocks the cold upstream matter can dissi- 
pate a significant fraction of kinetic energy, which is partially trans- 
ferred to a population of nonthermal relativistic electrons/positrons. 
Those particles emit synchrotron radiation and the most energetic 
of them are expected to cool rapidly enough, that the emission 
source is tightly localized in the shock vicinity. For a stationary 
relativistic shock, although the emitting elements can be moving 
with large Lorentz factor T = (1 — /3 2 ) -1 / 2 , (3 = v/c, the source 
position does not change in time. The relativistic enhancement of 
the intrinsically isotropic radiation is of the factor T> A /V (Sikora et 
al. 1997), where 



T (1-/3 cos £ fee ) 



(1) 



is the Doppler factor and £ fee is the angle between the direction 
of the element's motion and the observer direction in the external 
frame. Considering only the bolometric luminosity, we neglect the 
spectral index here. Thus, the Stokes parameters in the observer 
frame are: 



T) 3 



vQ v = II cos(2xe) vl v 
vTJ v = IT sin(2xfi) vl v ; 



(2) 

(3) 
(4) 



where II is the degree of linear polarization and xe is the electric 
vector positional angle (EVPA). Like in most studies of polariza- 
tion from astrophysical jets, we measure xe from the projected 
direction of the jet axis. 

In the following, we present our method of calculating II and 
Xe for emission from oblique shocks with chaotic magnetic fields. 
We choose a cartesian coordinate system (x, y, z), in which the jet 
axis is aligned with the z-axis and the observer is pointed by a unit 
vector k contained in the :rz-plane. The inclination of the observer 
direction to the jet axis is 9 obs . The shock front is inclined to the jet 
axis at angle r\ and its normal vector n makes an azimuthal angle 
(f> with the rrz-plane. A downstream fluid element is characterized 
with velocity j3 = v/c and the inclination of its velocity direction 
e to the jet axis 8 e . Thus, we have (see Fig. 1) 

k = [s'm9 obs ,0,cos9 obs ], (5) 

e = [sin 9 e cos 4>, sin 8 e sin (j>, cos 9 e ] , (6) 

n — [cos T) cos (j>, cos r\ sin 4>, — sin rj\ , (7) 

cos £fc e = k ■ e = sin 6 obs sin 9 e cos (j> + cos 6 obs cos 9 e . (8) 

We perform a Lorentz transformation into the fluid rest-frame 
O' . The observer position vector transforms like 

fc'=X>{fc + [(r-l)fc-e-rV?]e}. (9) 

The inclination of the shock surface transforms like 

tan(r7'-0 e )=rtan(77-0 e ). (10) 

Using 77', we find the normal vector to the shock surface: 




Figure 1. The coordinate system used for calculating the polarization de- 
gree and the electric vector positional angle (EVPA). The jet direction is 
along the z-axis. The observer is located in the :rz-plane, inclined to the 
jet direction by 8 obs . The black contour represents a shock surface element 
located at the positional angle <j>. Vector n, normal to the shock element, is 
inclined to the local radial coordinate by i). The downstream velocity field 
direction e makes an angle 6 e with the jet axis. Since we assume axisym- 
metric jet, vectors n and e are aligned with the <j> = const plane. 



n — [cost) cos 0, cos 77 sini 



- sin 77 



(11) 



Now we adopt a formula from Hughes et al. (1985) for the degree 
of polarization: 



n 



q + l (!-^) [l-(fc'-"') 2 ] 

a + 5/3 X 2-(1-k 2 ) [l-(fc'-n') 2 ] 



(12) 



where 1/n is the shock compression ratio and a is the spectral 
index of optically thin synchrotron radiation: F v oc p~ a . For the 
rest of this paper we will use a = 0.5, for which the maximum 
value of the polarization degree is Umax ~ 0.7. 

The magnetic vector of the polarized electromagnetic wave 
in O' is both normal to the propagation direction (observer) and 
tangent to the shock surface, thus 



B' oc ffe' x n') 



(13) 



We show in Appendix A, that the polarization angle is invariant 
in the Lorentz transformation, when measured with respect to the 
transformation vector, in this case e. First, we introduce an orthog- 
onal basis («', w') in the plane of the sky, for example 



v = [k' z ,0, -k' x ] , 



(14) 



© 0000 RAS, MNRAS 000, 000-000 



Polarization from reconfinement shocks 3 



w — [ k x ky, 1 hy , kyh^ . (15) 

Then, the positional angles of the wave magnetic field B' and the 
velocity direction e' = e are: 



tanxs 



B 'w> 

B' , 

v' 

tanx e = • 



(k' x n')-w' 
(k> x n>) ■ v' ' 



(16) 
(17) 



We go back to the O frame and find a basis, independent of that in 
the O' frame, in which the positional angle of the jet axis is 0: 



v = [-cos6» oi)s ,0, sm6 obs ], 
w = [0,-1,0]. 

The positional angle of the velocity vector is 

e w sin 9 e sin d> 
tanx e = — = -r-7 7. 7 Q ■ a — • 

e v sin 6 e cos 9 b s cos qb — cos 9 e sin 6J f, s 

Finally, the EVPA is given by 
2' 



XE = XB - Xe + Xe ~ 



(18) 
(19) 



(20) 



(21) 



For axisymmetric flows, the Stokes parameter U v averaged 
over the azimuthal angle <f>, for constant z coordinate, vanishes. 
The averaged EVPA can be either (parallel polarization, (Qv)^ > 
0) or 90° (perpendicular polarization, (Q„) ^ < 0). Thus we will 
use the quantity (Qv/h)^ as a polarization degree that contains 
information on the electric vector orientation. 



3 CONICAL SHOCKS 

CC90 studied conical shocks as a mean of deflecting a relativistic 
upstream flow parallel to the jet axis. They kept the angle between 
the shock surface and the jet axis r/ cc < tt/2. Thus, the matter 
crossed the shock from the outer side of the shock front (which 
is where the radial coordinate r is larger than the shock surface 
radius for given z-coordinate) to the inner side. If rf c > ty/2, 
the situation is opposite: the upstream side is the inner side and 
the downstream side is the outer side. It has been discussed in a 
later paper (Cawthorne 2006), that two conical shocks with rf( c + 
r/2 C = 7r and equal upstream velocities produces exactly the same 
Stokes parameters integrated over the azimuthal angle <p. This type 
of symmetry is also evident in our procedure described in §2. It 
can be easily shown, that transformation T = {r) — * —r\, 9 e — ► 
— 9 e , <j> ~ * 4> + 1"} preserves V, II and \e, so it also preserves the 
Stokes parameters. 

We calculate the polarization degree from conical shocks, 
when the upstream flow is not parallel, but diverges, making an 
angle 6 U with the jet axis. A model from CC90 with a shock incli- 
nation angle r\ cc and upstream velocity /3^ c is used for calculat- 
ing the downstream velocity value (3% , the downstream velocity 
deflection angle ^ cc and the shock compression ratio k cc from 
their Eqs. (1), (5) and (2), respectively. For the procedure described 
in §2 we take: jj = U - rf c , e = 9 U - £ cc , (3 = 0% C and 
k = k cc (see Fig. 2). For 9 U = this is equivalent to the CC90 
solution with <f> shifted by it. 

In Fig. 3 we present the averaged polarization degrees, as the 
functions of observer inclination o bs, for a series of models with 
cc = 10 and 9 U ranging from (CC90 solution) to 
i. e. a cylindrical shock surface). We find that for 



v cc = w ° r 
10° (r) = 0. 

observers with 9 b s < 25°, {Qv/Iv)^ strongly decreases with 
increasing 9 U , changing sign from positive (parallel polarization) 




Figure 2. Cross-sections through the conical shock: a) with parallel up- 
stream flow (like in CC90), b) with divergent upstream flow, where the 
shock surface (indicated by the thick gray dashed line), the upstream and 
downstream velocity fields (solid black lines) are all rotated by an angle 9 U 
with respect to the jet axis z. Thus, keeping the same value of upstream ve- 
locity (3u C , we obtain the same values for the downstream velocity f3j C 
and the compression ratio k cc . rj cc is the shock inclination to the up- 
stream velocity field and is the deflection angle between upstream and 
downstream velocity fields. The shock inclination r) is negative in both pan- 
els a) and b), and downstream velocity field inclination 8 e is negative in 
panel a), while positive in panel b), assuming 9 U > £p c ' . As discussed in 
the text, the structure shown in panel a) was obtained from a CC90 solution 
via transformation T. 
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Figure 3. Polarization degree from the conical shock averaged over the az- 
imuthal angle <f>, as seen by different observers, calculated for models de- 
scribed in the text and in Fig. 2, based on the CC90 solution with the shock 
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Plotted are the results for 8 U = (dashed line), 9 U € {2°, 4°, 6°, 8°} 
(gray lines) and 6 U = 10° (thick solid line). 
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Figure 5. Structure of the reconfinement shock for static external medium. 
The jet is launched at z = with half-opening angle Bj. The jet matter 
(upstream) crosses the shock front described by r s (z) and gets slightly de- 
flected in the downstream region. The shock reaches its maximum width 
r m at z = z m and ends in the recollimation point at z = z r ~ 2z m 
with the half-closing angle r ~ ©j. r c (z) denotes the contact disconti- 
nuity between the shocked jet matter and the external medium. This figure 
is taken from NS09. 

to negative (perpendicular polarization). For the cylindrical shock 
case a perpendicular polarization with degrees exceeding 20% can 
be easily obtained. For observers with o b s > 25°, (Qv/h)^ 
slightly increases with increasing 8 U , nevertheless the polarization 
degrees are very low. 

We have found, that further increase of 6 U in the series does 
not result in deeper perpendicular polarization. Thus, we confine 
ourselves to study the cylindrical shock case with varying 9 U . The 
results are shown in Fig. 4, which is analogous to Fig. 3b from 
CC90. We find, that the polarization is dominantly perpendicular, 
with the polarization degree not exceeding 28%. The maximum 
polarization degree is very similar in all models with 9 U > 20°, it 
is approximately obtained for 8 obs ~ 6 U . Only for small 8 U some 
observers would see a parallel polarization of very low degrees. 



4 RECONFINEMENT SHOCKS 

Following the idea that conical shocks can be approximated with 
a sequence of oblique shocks rotated around the jet axis (Lind & 
Blandford 1985), axisymmetric reconfinement shocks can be con- 
sidered as a sequence of conical shocks with different inclination 
angle 77. We use a semi-analytical model for the structure of recon- 
finement shocks from NS09 to provide downstream flow parame- 
ters for the polarization procedure described in §2. Basic param- 
eters of the reconfinement shock are presented in Fig. 5. Specif- 



Figure 6. Synthetic emission maps of the reconfinement shock with jet 
Lorentz factor Tj = 10 and half-opening angle ©j = 5° , as seen by differ- 
ent observers labelled with their inclination angle o b s . The jet propagates 
horizontally to the right direction: from the source marked by the cross to 
the recollimation point marked by the circle. The linear scale and pixel res- 
olution is the same in all cases. The gray shade is proportional to the total 
observed intensity /„, smoothed by a Gaussian distribution with standard 
deviation equal to the pixel diagonal size. The shading is normalized to the 
maximum pixel intensity calculated separately for each observer. Positional 
angle of the electric vector (EVPA) is indicated with bars in the center of 
each pixel. The bar lenghts are proportional to the polarization degree. A 
bar of length equal to the pixel size would correspond to a 100% polarized 
emission. The pixels of the lowest intensity were depolarized by adding a 
background of intensity equal to 3% of the maximum pixel intensity. 

ically, we take 'Model 2' described in that paper, which includes 
the transverse pressure gradient across the shocked matter zone, 
for the case of ultra-relativistic equation of state for the post-shock 
(downstream) matter (7s = 4/3). It is characterized by: uniform 
external pressure p e , bulk Lorentz factor of the jet Tj, half-opening 
angle of the jet Oj and total jet power Lj. As shown in NS09, the 
parameters p e and Lj have no influence on the shock shape, but 
they determine the size of the structure, which is represented by 
the reconfinement position z — z r . We assume that the luminosity 
of the emitted synchrotron radiation is proportional to the kinetic 
energy flux dissipated at the given shock surface element. 



4.1 Polarization maps 

In Fig. 6 we show synthetic maps of synchrotron emission from 
spatially resolved reconfinement shock, a Model 2 solution with 
Fj — 10 and Qj = 5°, as seen by observers located inside 
(dobs = ©j/2), at (8 bs = Oj) and outside (0 o bs = 29j) the jet 
opening cone. We find that they would all see an edge-brightened 
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Figure 7. Total observed emission intensity distribution along the jet di- 
rection, as seen by different observers labelled with their inclination angle 
9 bs ■ The profiles were obtained by integrating the emission maps similar to 
those shown in Fig. 6 across the coordinate perpendicular to the projected 
jet axis. They are normalized separately to their maximum value. The (- 
coordinate is defined in the text: £ = corresponds to the jet source and 
£ = 1 to the recollimation point. All models were calculated for jet Lorentz 
factor Tj = 10 and half-opening angle Qj = 5°. 

jet. There are two regions of enhanced brightness. One is located 
close to the jet source (marked with crosses), although it is much 
weaker for 9 a bs — 2Qj. The other is close to the recollimation 
point (marked with circles). These regions are characterized with 
the strongest Doppler boost, as for some azimuthal angles the shock 
surface is approximately tangent to the line of sight. 

The polarization degree maps are also edge-brightened, the 
polarization degree is especially small close to the jet axis. For 
small observer inclination (9 bs = Qj/2), when the jet appears 
only slightly elongated, the polarization vectors point approxi- 
mately radially outwards the midpoint between the jet source and 
the recollimation point. In the vicinity of the jet ends the polariza- 
tion is parallel (to the jet axis), while between them it is perpendic- 
ular. For observer located at the jet opening cone (9 obs = Qj) the 
polarization degree becomes lower, where polarization is parallel, 
while is high, where polarization is perpendicular. For large ob- 
server inclination (9 obs — 2Qj), the polarization is perpendicular 
everywhere in the structure. 

4.2 Longitudinal profiles 

When angular resolution is too low to resolve the transverse struc- 
ture of the jet, one may still obtain a longitudal brightness and po- 
larization profile. This is equivalent to integrating 2-dimensional 
maps like those shown in Fig. 6 across the coordinate orthogonal 
to the jet projected axis. To compare such profiles obtained for dif- 
ferently oriented observers, we introduce a coordinate £, which is 
related to the coordinates used in Fig. 1 by 

C = . (22) 

For points on the jet axis x = 0, thus £ = z/z r . For the recollima- 
tion point C = 1- 

In Fig. 7 we show the profiles of the total observed intensity 
for the same structure as in §4. 1 . Each profile has a clear maximum, 
which is close to the recollimation point for small 6 obs (£ ~ 0.85 
for 9 bs = Qj/2) and shifts towards the jet midpoint as o b s in- 



Figure 8. Polarization degree profiles along the jet direction, as seen by 
different observers labelled with their inclination angle o b s - They were 
obtained from the Stokes parameters Q v and I v integrated in the way de- 
scribed in Fig. 7. All models were calculated for jet Lorentz factor Tj = 10 
and half-opening angle Qj = 5°. 



creases (£ ~ 0.6 for o b s = 30°). A secondary maximum ap- 
pears close to the jet source (£ ~ 0.2) for 6 b s smaller (and ap- 
parently also for slightly larger) than Qj . Thus the profile changes 
from double-peaked to approximately parabolic, as o b s increases. 
This is consistent with total intensity maps from Fig. 6. This profile 
morphology change reflects the importance of relativistic boosting 
of the radiation. As was shown in Fig. 9 of NS09, most of the ki- 
netic energy dissipation takes place for C, ~ 0.5. The luminosity 
profiles for large o b 3 are consistent with this, but the profiles for 
dobs < Qj show effects of strong Doppler boosting from the jet 
portions characterized with \r/\ ~ 9 bs, thus producing peaks that 
do not correspond to the position of maximum kinetic energy dis- 
sipation. 

In Fig. 8 we show the polarization degree profiles for the same 
structure. Observers with o bs < Qj would see parallel polariza- 
tion from regions with ( < 0.1 and £ > 0.9 and perpendicular 
otherwise. The maximum parallel polarization degree is obtained 
outside both projected ends of the structure (C < and ( > 1) and 
is higher for 6 ob s = Qj/2 (~ 40%) than for 6 ob s = Qj (~ 20%). 
The maximum perpendicular polarization degree (~ 40%) corre- 
sponds to ( G [0.4; 0.6]. 

For 9obs > Qj the observers see perpendicular polarization 
for all £. The polarization is only significant between the ends of 
the structure, namely for < £ < 1. The maximum polarization 
degree is decreasing with increasing 9 bs, changing from ~ —35% 
for 9 ob s = 2Qj to ~ -10% for 6 obs = 30°. The polarization 
degree profile changes from approximately symmetric with respect 
to £ = 0.5 for 9 bs = 20j to that with the polarization degree 
increasing linearly with £ for 9 a bs = 30° . 

Note that all of the observers see perpendicular polarization 
from the region with £ G [0.1; 0.9]. In this region they see dom- 
inantly the portions of the shock with small inclinations |rj| < 
Qj/2. Thus, the perpendicular polarization seen there is consis- 
tent with the cylindrical shock solutions discussed in §3. As one 
can see in Fig. 7, this is also the region, from which the bulk of 
the total emission originates. Thus, we expect that this region dom- 
inates the radiative output from the reconfinement shock and so the 
averaged polarization from the whole structure is perpendicular for 
most observers. 
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Figure 9. Spatially averaged polarization degrees from the reconfinement 
shock, as seen by different observers labelled with their inclination angle 
9 obs . Models were calculated for jet Lorentz factor Tj = 10 and different 
jet half-opening angles Oj . 
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Figure 10. Spatially averaged polarization degrees from the reconfinement 
shock, as seen by different observers labelled with their inclination angle 
9 obs . Models were calculated for jet half-opening angle Oj = 5° (marked 
with the gray vertical line) and different jet bulk Lorentz factors Fj . 

4.3 Spatially averaged polarization 

In Fig. 9 we show the polarization degree averaged over the whole 
spatial extent of the reconfinement shock for four models with 
Tj — 10 and different half-opening angles Qj. The polarization 
is dominantly perpendicular, with maximum degree observed for 
dobs > Oj. The maximum polarization degree value is decreas- 
ing with Qj, ranging from ~ 27% for Qj = 2° to ~ 13% for 
Qj — 20°. For Qj > TJ 1 observers closely aligned with the jet 
axis will see some parallel polarization, with maximum polariza- 
tion degree ~ 7% for Qj = 20°. For observers with large 6 ooa the 
polarization degree increases with Qj . 

In Fig. 10 we show the spatially averaged polarization de- 
gree from the reconfinement shock for four models with Qj = 5° 
and different bulk Lorentz factors Tj. Again, the polarization is 
mostly perpendicular, with the exception of closely aligned ob- 
servers, when Tj > QJ 1 (with maximum parallel polarization de- 
gree ~ 7% for Tj = 40). The position of maximum polarization 
degree approaches Qj with increasing Tj and the maximum polar- 
ization degree value decreases, ranging from ~ 26% for Tj = 5 to 



~ 14% for Tj — 40. For observers with large 6 ooa the polarization 
degree significantly decreases with r 3 . 



5 DISCUSSION & CONCLUSIONS 

We have shown that the degree of perpendicular polarization from 
relativistic conical shocks with chaotic magnetic fields can be large, 
if one considers divergent upstream velocity field. In the extreme 
case, when the shock surface is parallel to the jet axis (a cylindri- 
cal jet) the degrees of perpendicular polarization may reach 28%. 
Thus, the conical shocks can account for the observed polariza- 
tion bimodality of parsec-scale jets on mm wavelengths (Nartallo 
et al. 1998; Jorstad et al. 2007), without introducing parallel mag- 
netic fields, as proposed in Cawthorne (2006). It is important to 
notice, that one can only consider 'local' solutions of this kind, i. e. 
axisymmetric sections that are far enough from the axis, to avoid 
problems on the r — boundary. 

Our results on polarization from the conical shocks have im- 
portant implications for the reconfinement shocks, that are mod- 
elled in terms of an interaction between the spherically divergent 
relativistic jet (upstream) flow and the static external medium. Ev- 
ery particular section of approximately paraboloidal shock surface 
can be obtained from a CC90 solution via simple rotation described 
in Fig. 2. The main difference between our and their models of the 
shock configuration is that we assume a cold upstream matter, mak- 
ing possible the solutions for arbitrarily small inclinations of the 
shock front, with respect to the upstream velocity field. 

We have found, that the emission from axisymmetric recon- 
finement shocks with chaotic magnetic fields is clearly dominated 
by perpendicular polarization. Both total intensity and polarization 
degrees are higher on the edges of the structure and the wave elec- 
tric vectors are perpendicular to the outline of the shock surface. 
The polarization maps are axially symmetric, in accordance with 
the absence of a large-scale helical magnetic field component (Lyu- 
tikov et al. 2005). The total intensity maps are strongly affected by 
Doppler boosting for observers located closely to the jet opening 
cone. One can see two 'knots' for o b s £ {©j/2, Qj} close to 
both of the shock ending points, but only one knot for 9 obs = 2Qj 
close to the recollimation point. On the longitudinal profiles (Fig. 
7) we observe this knot to shift towards the shock midpoint and 
loose its brightness contrast. The polarization of the knots can be 
estimated, by summing the Stokes parameters from the 3x3 pixel 
groups taken from Fig. 6 and centered on the brightest pixel of the 
knot. It is parallel with degrees about 30% for both knots seen for 
dots = ©j/2. The knots seen for 9 obs = Oj are effectively depo- 
larized, the one close to the source of the jet having 4% net paral- 
lel polarization, while the other one showing 1% net perpendicular 
polarization. Finally, the knot seen for o bs = 20 j shows perpen- 
dicular polarization with degree of 23%. Thus we predict that the 
polarization of the knots associated with the reconfinement shocks 
may be both parallel or perpendicular, depending on the observer 
position with respect to the jet opening cone. 

The longitudinal profiles of the polarization degree (Fig. 8) 
show that observers closely aligned with the jet axis may see po- 
larization degrees as high as 40%, both parallel and perpendicular. 
This is higher than the values obtained for conical shocks. It can be 
explained as a projection effect, since points of equal £ coordinate 
form a plane oblique to the jet axis. Integrating emission from such 
planes includes shock surface portions with different inclination to 
the jet axis. If we instead integrated emission from the points of 
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equal z coordinate, we would obtain polarization degrees consis- 
tent with the solutions for conical shocks. 

On the longitudinal profiles for closely aligned observers, the 
transition between the parallel and perpendicular polarization is 
approximately coincident with the total intensity profile maxima. 
Note, that for o b 3 — Oj/2 the total intensity profile maxima are 
slightly shifted towards the jet midpoint from the positions of the 
highly parallel-polarized knots. 

The spatially averaged polarization degrees also prove 
the dominance of perpendicular polarization. The models with 
TjQj < 1 produce perpendicular polarization for all observers, 
while the models with FjOj > 1 can produce parallel polariza- 
tion only for the observers located inside the jet opening cone. This 
results from very strong Doppler boosting from the shock surface 
regions tangential to the line of sight. The transition between polar- 
ization degree maximum and minimum always takes place around 
dots = ©j and its sharpness —d({Q v / 7„) ^ z ) / 'd8 obs increases 
with Tj. Thus, extremely relativistic reconfinement shocks could 
produce the polarization angle swings with moderate change of jet 
orientation. 
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APPENDIX A: LORENTZ INVARIANCE OF THE 
POLARIZATION ANGLE 

It is well known that the polarization angle of linearly polarized ra- 
diation is invariant in Lorentz boost, if this angle is measured from 
the plane containing both the wave propagation direction and the 
boost velocity vector (Cocke & Holm 1972). Straightforward for- 
mulae for the value of polarization angle can be derived, when one 
chooses to set up his coordinate system aligned with this preferred 
plane (Blandford & Konigl 1979; Bjornsson 1982). However it is 
sometimes more convenient to use a fixed coordinate system, in 
which the wave propagation direction and the boost direction are 
completely arbitrary. One has also a freedom of choosing a basis in 
the plane of the sky, in which polarization angle is associated with 
the azimuthal angle of the wave electric vector. The formulae for 
the synchrotron polarization angle for arbitrary basis and magnetic 
field direction have been given by Lyutikov et al. (2003) (see their 
Appendix C). If the plane-of-the-sky basis is chosen to be somehow 
aligned with the fixed coordinate system, the polarization angle has 
to be transformed in non-trivial way. But in this transformation it 
this also the basis that rotates, and so the projection of the velocity 
vector onto the basis is different. Below we show in most general 
approach, that the angle between the polarization vector and the 
projection of the transformation velocity vector onto the plane of 
the sky is Lorentz invariant. 

Consider a frame O, in which a linearly polarized wave prop- 
agates towards an observer pointed by unit vector k. Frame O' is 
moving with the velocity (3 = j3n, where \n\ = 1. The polariza- 
tion vector positional angle x of the wave is measured in the plane 
of the sky, begining from n projected onto the plane of the sky. We 



construct an orthogonal basis (v, w) in the plane of the sky, that 
fulfills the conditions n • v > and n ■ w = 0: 



n — (n ■ k)k 
v = = — i 
\A-(«-fc) 2 

k x n 

w = k x v = — 

The electric vector of the wave E 
mined by the polarization angle: 

e = cos x v + sin x w. 



(Al) 
(A2) 

Ee has its direction deter- 
(A3) 



The magnetic vector of the wave is B = Bb, where B — E and 
b = k x e. The Lorentz transformation of electric field is given by 
(Rybicki & Lightman 1979) 



E'\\ = E h 

E' ± = r(E ± +(3xB). 

Hence the electric vector in O' is 

E' = TE - (r - l)(n ■ E)n + T/3n x B. 

One can show, that 

E 



E' = T(1- f3n-k)E = 



V 



The field direction can be written as 

, _ cos x I + sin x k x n 

where 

l=[l-2>(r-l) (l-(n-fc) 2 )] 



n • k — j3 



1-Bn-k 



(A4) 
(A5) 

(A6) 
(A7) 
(A8) 

k. (A9) 



Note, that I ■ (k x n) = and \l\ = y/l - (n ■ k) 2 , so in Eq. 
(A8) we have effectively decomposed the transformed electric field 
direction in a new orthogonal basis. Thus, we conclude that 

tanx'=tanx, (A10) 
i. e. the polarization angles are the same in O and O' . 
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